Conservation Laws and Cosmological Perturbations in 

Curved Universes 



Jean-Philippe Uzan\ Nathalie Deruelle^'^ and Neil Turok^ 
^ Departement d'Astrophysique Relativiste et de Cosmologie 
UPR 176 du Centre National de la Recherche Scientifique 
Observatoire de Paris, 92195 Meudon, France 

^ Department of Applied Mathematics and Theoretical Physics 
University of Cambridge, Silver Street 
Cambridge CBS 9EW, England 

February 7, 2008 

Abstract 

When working in synchronous gauges, pseudo-tensor conservation laws are often 
used to set the initial conditions for cosmological scalar perturbations, when those 
are generated by topological defects which suddenly appear in an up to then per- 
fectly homogeneous and isotropic universe. However those conservation laws are 
restricted to spatially flat {K = 0) Friedmann-Lemaitre spacetimes. In this paper, 
we first show that in fact they implement a matching condition between the pre- 
and post- transition eras and, in doing so, we are able to generalize them and set 
the initial conditions for all K. Finally, in the long wavelength limit, we encode 
them into a vector conservation law having a well-defined geometrical meaning. 

Pacs : 98.80.Hw, 98.70.Vc, 04.20.Cv 
Ref. PRD : DD6306 

1 Introduction 

Topological defects act in cosmology as "active" sources for the metric and matter per- 
turbations in that the linear differential equations for the evolution of these perturbations 
are inhomogeneous (in contrast with what happens in inflationary scenarios). Thus the 
general solution of the evolution equation for each perturbation mode is the sum of a 
particular solution of the complete equation and the general solution of the homogeneous 
equation. It depends on integration constants that have to be fixed by the physics which 
gave rise to the defects. 

One has thus to take into account the fact that the defects appear at a phase transition 
in a previously strictly homogeneous and isotropic universe. For scales much larger than 
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the Hubble radius at that time, which are the scales of interest today, the phase transition 
appears as instantaneous (in the following we shall use indifferently the words "scales" or 
"wavelengths"). A way to fix the initial conditions is therefore to find linear combinations 
of the perturbation variables which are constant in time and force them to vanish since 
they were strictly zero before the phase transition. It has been shown [||] that, in the 
newtonian gauge, the initial conditions are completely fixed by these matching conditions. 

In synchronous gauges, for flat {K = 0) Friedmann-Lemaitre spacetimes, and for 
"scalar" perturbations (see [0] for the decomposition of perturbations into scalar, vector 
and tensor modes), the initial conditions are usually fixed, at least partially (see e.g. 
@|) 0) 111) d]))) by means of a "conserved pseudo-tensor" r^,^ more precisely by the 
constraint "roo = 0" (the equation of conservation doToo = dkTok — on large scales then 
ensures that this constraint is conserved during the evolution [^). 

In this article, we first compare and contrast the implementation of the initial condi- 
tions for long wavelength scalar perturbations in both gauges (& 2 and & 3). To do so 
we consider a toy model, namely coherent defects in a radiation dominated, spatially fiat, 
universe. We stress some of the Charybdes and Scyllas of studying its perturbations in 
synchronous gauges, and show how the "pseudo-tensor" r^i, defined in partially fixes 
the initial conditions. We then interpret the condition "roo = 0" as a matching condition 
between the pre- and post- transition eras. Finally, we show how the newtonian gauge 
simplifies the resolution of the evolution equations. In & 4, because we had related it to 
a matching condition in the fiat case, we can construct a generalized "conserved pseudo- 
tensor" in synchronous gauges for K Friedmann-Lemaitre spacetimes. Finally, in & 
5, we formalize the definition of this "pseudo-tensor" by means of a well defined vector 
conservation law [j^ and give it a precise geometrical meaning. 

2 Scalar perturbations in synchronous gauges 

The line element of a perturbed Friedmann-Lemaitre spacetime reads, when the pertur- 
bations are scalar and in synchronous gauges : 

ds"^ = 0^(1]) 

where t] is conformal time, a{r]) the background scale factor ; 'jij is the three-metric of the 
constant time maximally symmetric surfaces in the coordinate system (i = 1, 2, 3) ; Di 
is the associated covariant derivative, A = DiD^ (latin indices are raised by means of the 
inverse metric 7*-') ; finally h and E are two "small" functions of space and time. 

Infinitesimal coordinate transformations which preserve the synchronicity of the gauge 
can be performed. They are defined by : 

r]^r] + T , x''^x'' + d''L (2) 

where T and L are two first order functions of rj and such that 

T' + nT = and L' - T = ^ T = — and L = d + eo / — (3) 

a J a 



-drf + 



;i + 3^)7ii + 2(A/^i - ^liA)E \ dx'dx^ 
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where a prime denotes a derivative with respect to conformal time, where Ti = a' /a and 
where cq and ei are functions of space only. The two scalar metric perturbations h and 
E then transform as 

h^h + QnT + 2AL , E^E + L (4) 

with T and L given by (3). 

The energy-momentum tensor of the matter content of this perturbed universe can be 
written as : 



T. 



(5) 

Tfj,u (/U, ^ = 0, 1,2,3) is the energy-momentum tensor of the homogeneous and isotropic 
background of radiation and dust ; ST^^i, is its perturbation : in synchronous gauges, 
its scalar components can be expressed in terms of two scalar quantities, 6 = Sp/p, the 
density contrast, and v, the velocity perturbation, as (see e.g. IQ) : 



ST, 



00 



pa^S, STqi = — pa^(l + uj)DiV, ST^ = Pa^ 



^ + r + ^5 ) + 2DiDjE 



(6) 

p and P are the density and pressure of the background fluid, u = P/ p, cj. = P'/ p' and 
PT = Sp- cl5p. 

Finally 0^,^ is the energy-momentum tensor of the topological defects. We suppose 
that it is a small perturbation which does not contribute to the background. We decom- 
pose its scalar components as : 



Boo = p", Ooi = -Div' 



7., P 



The four source functions p^,P*,f'*,n* will be discussed later. 

The energy and momentum conservation equations for B^j, then read : 

p'' + n{p' + 3P') + Av' = 

v'' + 2W + P' + -Aff = 0. 

3 

The conservation equations for the fluid are, at linear order (see e.g. 0) : 



(7) 



(9) 



5' + -{l + uj)h' 



[l + uo)Av 



(10) 



-5. 



v' + n{i-?.cl)v 

As for the linearized Einstein equations they can be cast under the form (see e.g. 
0): 



Hh' = K{a^p6 + p') 

a'^p{l + uj)Av + Av 
h" + Hh' = -K aV(l + 3c^)5 + {p' + 3P' 



^A + ir) hr' -Kh' 



K 



(11) 

I or 
(12) 

(13) 
(14) 
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{h - h^y + 2n{h - h-y - -Ah- = -2kau' (is) 

3 

where k = 8ttG, G being Newton's constant, and where the perturbation h~ is defined as 

h- = h-2AE. (16) 

Equations (|1-|15D are eight equations for the four unknowns h,h (or E), 6 and v : 
the source functions p*, P"^, v'^, H**, subject to the constraints (H-PD being known, equations 



( p!0| - |T5|) must give h,h~ (or E), 6 and v, two of these equations being redundant. The 
four Einstein equations for scalar perturbations ([T^^) include two constraints on the 
four variables so that the general solution depends on two physical degrees of freedom 
which must be fixed by the physics of the problem. 

To be complete, let us recall the Friedmann equations for the background 

Ka^p = 3{n'^ + k), Ka^p + 3P) = -GH'. (17) 

3 Synchronous gauges inside out 

3.1 A toy model : coherent defects in a flat, radiation domi- 
nated, universe 

To show, on a simple example, how the above system of equations can be solved, we 
consider here a radiation dominated universe with fiat, K = 0, spatial sections. Hence, 
u = = 1/3, T = and, from (0), 0(77) oc rj. Moreover we suppose that we are deep 
enough in the radiation era for the perturbations of interest today to be all larger than 
the comoving Hubble radius {l/Ti). Finally we assume that, in that large scale (or long 
wavelength) limit, the variables describing the defects take the form 

p« = r]-^/^Ai, P' = r/-i/2^2, V = -v^/^As, W = if^^A^ (18) 

where Ai are four functions (or, rather, four random variables) of space varying on scales 
larger than the Hubble radius. Provided that the defect network is coherent, uncorrelated 
on scales larger than the Hubble radius and scales with the background (the precise 
definition of these properties is given in e.g. |I|), it can be modelled by (|TB|), (see e.g. 
|jlO|). Now, the precise energy-momentum tensor of the network of topological defects 



depends on the kind of defects considered and can be obtained only by heavy numerical 
simulations (see e.g. 0). It turns out that, generically, the network is not perfectly 
coherent so that AH'*, and not H'* itself, varies on scales larger than the Hubble radius. 
However such incoherent sources can, in principle, be decomposed into a sum of coherent 



sources described by ( p^q ) ||T3[| . 



The conservation equations (§-0) first impose that, in the long wavelength limit (where 
we neglect the gradient terms) 

= -QA2 and A3 = Ia2. (19) 

5 



Then equation (|10]) and (|T^) can be rewritten, in the same approximation, using (|T8D and 
(pjp, as an evolution equation for 5, 

5" + -5' - ^5 = -2KA2r]-^'\ (20) 
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the general solution of which is 



6 = W + ^ + ^«:^2r/^/^ (21) 

where Aq and Cq are two constants of integration (that is two slowly varying functions 
of space). It is usual to call the Aq term the "growing" mode 0]. Then the continuity 
equation ([T0|) yields 

h = -h,r^' - 6^ - ^4^A,r^'^" " 9^o, (22) 
2 r]^ 7 

\l/o being a new function of integration. The Euler equation ([TT|) then gives v, 

Ai being another function of integration. 

At this stage of the resolution we are left with the problem of determining the sec- 
ond scalar perturbation of the metric (either h~ or E) with one of the three remaining 
equations (P^, (|T3|), (|T3p. For instance, the resolution of (^) yields 

E = -e^-eo\nv + ^Xov' -^v'- — + ^^ (7^ - 2^) r/^/^, (24) 
80 A rj 441 

ei and being yet two other integration functions. 

On the six integration functions that have been introduced, namely Aq, cq, ^E'o, Ai, ei, 
and \E'i, two of them (cq and Ci) can be eliminated by a gauge transformation (see eq. H). 
This means that the two remaining Einstein equations (|T2[) and (|I3|) are not identically 
satisfied, but, rather, should give two constraints on the remaining four integration con- 
stants, leaving out the two physical degrees of freedom. This task requires however some 
shrewdness since the gradient terms in ([l2|) and (|l^) cannot be simply disposed of. 

In another route was therefore pursued. First a "conserved pseudo-tensor" r^i, 
defined by [|lTl 

KToo = K {6Too + Qoo)-Hh', KTok = K{6Tok + Qok), t^Tki = K{6Tki + Qki) -^{h'j^i-h' 5ki) 

(25) 

was introduced. This quantity is indeed conserved, in that Tq^ = djTij and Tqq = djTQj 
provided Einstein's equations (8-17) are satisfied. It has then been widely used (see e.g. 
0) 0) 0) 0) to fix the initial conditions, that is to determine the value of some of the 
arbitrary functions of space previously introduced. The line of the argument is that since 
the defect network is uncorrected on scales larger than the Hubble radius, Tki must be 
white noise on large scales, that is that its components, like the source functions Ai, must 
be functions of space varying on scales larger than the Hubble radius. Now, from the 
conservation laws one gets that Tqq = dijTij. Hence, if Too and Tqq are continous through 
the phase transition, then Too is second order in the gradients, i.e. 

Too = (26) 
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and the expressions 



on large scales. Now, from the definition 
Qfj.u, we have 

3. 6A2 h' 
Too = ^0 — = oAo 



(0) for 6T^u and 
(27) 



(using the solution (|21| ) and ( [2^ ) for S and h). Therefore the integration function Aq 
is second order in the gradients and must be set equal to zero. Thus the use of the 
"conserved pseudo-tensor" (25) eliminates the "growing" mode of the density constrast 
6. 

As for the other integration functions \E'o, ^^i and Ai entering the solutions (|2^), (p3| ) 
and (24) for h, v and E (cq and ei just refiect the non-unicity of synchronous gauges) 
they are usually ^ chosen by assuming continuity through the phase transition, that is 
set equal to zero. 



3.2 Matching conditions in synchronous gauges 

The choices made above to fix the functions of integration (that is to impose the continuity 
of Too, h and v and E) must be justified. The question is therefore to determine which 
quantities are indeed continuous through the phase transition. This has been studied 
in detail in For perturbations varying on scales larger than the Hubble radius the 
phase transition which is at their origin looks instantaneous and the matching conditions 
between the pre- and post- transition eras are the standard ones : the induced three- 
metric on, and the extrinsic curvature of the surface of transition, which is taken to be a 
surface of constant density, must be continuous. 

In synchronous gauges, these matching conditions translate as (we give them for 
all K) 



h- + 
HE' + 



Hh-' - 1 - - 



3in^ + K){l + iu] 

1 

9(n^ + K){l + u] 
K 



K0?p5 + Kp' 







3(H2 + ir)(l + cu) 



'p6 + np'^ 



(28) 
(29) 
(30) 



where [F]± = lim^_,Q+ {F{riPT + e) — F{riPT — e)), ?7pr being the conformal time at which 
the transition occurs. Two useful linear combinations will turn out to be 



'h- - GHE' 



J ± 







0. 



(31) 
(32) 



3.3 Pseudo-tensor and matching conditions 

Let us now study the implications of these general matching conditions when applied to 
our toy model. 
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First, equation (0), when particularized to a radiation background, large scale fluc- 
tuations and flat spatial sections, is nothing else than the condition of continuity of Too- 
This justifies the assumption (p6D and hence the choice Aq = (see (^)). 

In order now to assess the validity of the other choices made (continuity of h, v and 
E), we shall turn to the newtonian gauge. Indeed, (|3^) implies, in our toy model, that 



h- 



-E' 



(33) 



J ± 



and thus that, see (p!6|) 

h- = OiE), h = Oih-). (34) 

A natural set of variables to solve Einstein's equations (8-15) is then {h^,E) or {h,E), 
but not {{h — h~),h). Now, it also makes sense to find a combination of the perturbations 
variables such that the set of equations ( [Tol - p^ ) reduces to a master equation and to 



five algebraic equations, two of which being truly redundant, and not, as is the case in 
synchronous gauges, constraints on integration constants. Such a combination is known 
and amounts to working in the newtonian gauge H]. Indeed, defining 



h' 



HE', $ = -E" - HE', 



6 + 3n{l + uj)E', 



V + E\ (35) 



(with H = 1/r] and = 1/3 for our toy model) it is easy to show that (p!0| - p!5D reduce to 
(see e.g. ^) : 

3 



3 



u"' + $ + 
^ _ $ 

V V 
1 



0. 



K 

2 



V 



7]^ 



4 1 / 1 n'' 1 1 

^" + - -A^ = K -Aff + — + -P' - -p' 



(36) 

(37) 
(38) 

(39) 
(40) 
(41) 



The sources being known, the integration of (^) gives \E' in function of two integration 
constants. Then $, f", and 5" follow algebraically from (^), (^OD, and (|39D , that is 
without introducing new constants. More precisely the solution of this system is (see pi) 



^1 



9 



^0 + ^ + ^/^^/3/2(2^2-7^4), 



9 



(42) 



2^o + 4^ + ^/€r7=^/'(A2 + A4) 



+ ^ - l^r]^/\4A,-5A,). (43) 
2 r]^ 45 



rj- 9 

Only two integration "constants", \E'o and \E'i, are introduced (and not six as is the case 
when working in synchronous gauges), so that the conservation equations ( |5BD and ( |57| ) 
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are identically satisfied. Note that in the newtonian gauge, the "growing" mode is the 
term proportional to \l/o (and is not the same as the synchronous growing mode which is 
proportional to Aq, see (PD). 

Now the relations (|35| ) can be inverted to compute h,h~,S and v. Since one has to 
solve two differential equations to do so, two functions of space appear, which can be 
identified with cq and ei, and the comparison of (^) and ( ^5D to ( pij - piD tells us that 



Ai = and Aq = 0, (44) 

at lowest order in the gradients, which justifies the assumption of continuity of v, and, 
again. Too (tqo = for large scales being nothing more than a rewriting of (p9|)). 

Finally, the two constants \l/o and correspond to the two physical degrees of freedom 
and are determined by the matching conditions. They turn out not to be zero, but the 
terms proportional to them in (42), (43) become nevertheless quickly negligible (see [|l]). 
Hence the assumptions that h and E be continuous are therefore asymptotically valid. 

To conclude, working in synchronous gauges introduces additional integration func- 
tions which render the resolution of the perturbation equations slightly cumbersome. A 
solution to get rid of these spurious functions is to go to the newtonian gauge where 
the perturbation variables are such that the set of equations (|IU[|T^) reduces to a unique 



master equation and algebraic equations, two of which are redundant if the others are 
satisfied. The matching conditions at the surface of transition then fixes the two physical 
functions and justify the use of Too to eliminate the synchronous "growing mode" (i.e. Aq), 
and the assumptions of continuity of v and (approximately) h and E. Now, no gauge is 
perfect : when the sources are incoherent, 11'' in the right-hand-side of (ffTf) is a non local 
term. 



4 Generalization to universes with hyperbolic sec- 
tions 

We now turn to universes with hyperbolic hypersurfaces [K = — 1). We will first work in 
the newtonian gauge and then generalize the concept of "pseudo-tensor" . 

When the spatial sections are hyperbolic, the curvature generates a new length scale 
on top of the comoving Hubble or horizon scale [Rh = ^fH), which is the comoving 
curvature scale {Rc = !)• As follows from ([T7|) , the curvature scale is larger than the 
horizon scale 

Rh < Rc- (45) 

We are interested in all scales L that were larger than the horizon scale at the phase 
transition. They fall in three categories : (I) those which entered the horizon between 
the phase transition and today (i.e. L < Rro), (H) those which have not yet entered the 
horizon but are subcurvature (i.e. Rhq < L < Rc), (HI) those which are supercurvature 
(i.e. L > Rc). 

4.1 Newtonian gauge 

As we have seen, the introduction of the newtonian gauge perturbations ( p5|) was useful 
to keep to a minimum the number of integration functions. The generalization to 7^ 
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spacetimes governed by a mixture of fluids of the master equation for \1/ (^Tj) is (see e.g. 
i) 

^" + 3(1 + cl)W + [2n' + (7^^ - K){1 + 3c^)]^ - cl^m = S (46) 
where 5 is a source involving the variables describing the topological defects, whether 



they satisfy ([T9| ) or not. Equation (|46|) is closed, in the sense that 5 is a known function 
of time, only in the case when the background is governed by a single fluid. In the general 
case of a multi-fluid system, one has to add a set of equations for the relative entropies 
which enter S (see [H). However, on super-horizon scales, all these entropies turn out 
to be constant so that the multi-fluid can be described as a single one and the equation 
is again closed. In this paragraph, we therefore restrict our attention to modes which 
have not yet entered the horizon that is such that (L > Rhq), i-e. from categories II and 
III. For those modes, it is easy to show that the term c^A\l/ is negligible compared to 
(7^^ — K)'$ in the regime when Ti^ is of order K (see Appendix 1). 
Equation (^) can then be rewritten under quite a compact form as 

where we have introduced instead of the new variable x defined by 

1 (^a'\ , . 



The general solution of (^) can be written formally as 



X = Xo + I dn ^^ ^, , (49) 



where Xo is an integration constant, that is a slowly varying function of space only. 
Inverting equation (|48|), we can therefore write the general solution of equation (|46| ) on 
superhorizon scales as 

^^ = X.-Xoj dr^a jdr^a jdv ^,_^,^j^ S (50) 

where xi is another integration function, xo and Xi are determined by the matching 
conditions as in the flat case (in which case they are proportional to and \E'i). 

(Note that in the absence of the source term 5, as is the case in inflationary models, 
X generalizes Lyth's variable C, to Friedmann-Lemaitre spacetimes with hyperbolic 
spatial sections. Its constancy on superhorizon scales (equation (^9])) can then be used to 
compute, by means of (0), the amplification of ^ in open inflationary models.) 

As for the other perturbations, they are given algebraically by 

2 / K 

{k{p^ -3W) -QH{^' + n^))}. (51) 



3(^2 + K) 



(which generalize equations (p8|) (Wi) (HOl) (^T])). The solution depends only on the two 



integration functions Xo and Xi? and the continuity and the Euler equations are identically 
satisfied. 
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4.2 Synchronous gauges 



Provided we have solved the equations of the background (]T7|), the set of equations |n|, 



141) can be solved as in the section |3]T|. However, the question is then to find a prescription 
to fix the coefficient of the synchronous "growing mode" of the density contrast (i.e. the 
equivalent of the constant Aq in (pTD). For that purpose, we will generalize Too- 

If we note that, when K = 0, tqo = 0, Too being given by (^), is in fact a way to write 



the matching condition (|3T| ) at lowest order, a straightforward generalization of Too is 

KToo = ^{K6Too + Keoo + Kh- -Hh'} (52) 
where 7 is the determinant of the 3-metric 'jij. Indeed, if we introduce 

KTok = {n5Tok - 2KduE'} , (53) 



it can easily be checked, using the Einstein equations (|T2], |T3|), that Toy is conserved, in 
the sense that 

doToQ = dkTok- (54) 
Now, the matching condition (|3TD imposes that initially 

Too = (55) 



at lowest order on superhorizon scales. Furthermore ( |5^ ) implies that Too remains zero at 
all times for superhorizon modes. Thus ( ^5]) is a way to eliminate at all times the spurious 
growing mode in the solution for the density contrast 6. 



5 The matching conditions in terms of a vector con- 
servation law 

In the previous section we have encoded part of the matching conditions into a quantity 
Too, the geometrical status of which is not defined. We show here that for long wavelength 
modes, it can be interpreted in terms of a conserved vector related to a symmetry of the 
background. For that purpose, we briefly describe the formalism we use and then apply 
it to our problem. 



5.1 General Formalism 

A general formalism was developped by Katz-Bicak-Lynden Bell |^ to define conserved 
quantities and conservation laws in an arbitrary spacetime {M-ig^i)- This generalizes the 



work by Bergmann who defined conservation laws with respect to a flat background |l5 . 
The formalism involves the introduction of a background {M.,g^jj) and a mapping, i.e. a 
way to identify points of }A and M.. 

The central idea is that one can construct a lagrangian for the gravitational field, 
quadratic in the covariant derivatives of the metrics and normalized so that it vanishes 
on the background (i.e. when M. = M), which is a true scalar density. From this 
quantity, one can therefore build vectors densities which are conserved in the sense 
that dj^ = 0. 
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The detailed construction of the vector densities and the demonstration of their 
properties can be found in [0. We will just sum up the relations which will be useful and 
fix the notations. The conserved vector field density is given by 

1 

— L 

2 



where is an arbitrary vector field. For any quantity A, A denotes y/—gA, g being the 
determinant of the metric g^,y, and A the value of A on the background. and are 
the Einstein tensor densities of the spacetime and of the background, l'^'^ is the difference 
of the two metric tensor densities, 

R^i, is the Ricci tensor of the background spacetime. Introducing 

where F^^ and F^^ are the Christoffel symbols of the spacetime and of the background, 
as well as and the two covariant derivatives associated with g^y and gfj,^, the 
expressions for t^, a^^P'^^ and Z^^^") are 



2Kt> = q'"^fA\A'' +A''A^ -2A^A^ 



+3^" (ALAj. - A^^Ai) - r (A^^A^^ - A^.AJ,) 5^, (57) 

(this term reduces to the Einstein pseudo-tensor density when the background is Minkowski 
spacetime in cartesian coordinates), 

2/s:<T^[^"] = - g^H''^^'^ A^^ - 2Z^[^^"1'^A^^ (58) 

(where [] means antisymmetrization), and 

4KZ^(r) = {z^g''' + g^'Z;-g'^-z)A^^^ 

+ {g^-Z - 2r^Zl) A;^, + r^dxZ + r (V^Zp. - 2VpZ^) , (59) 

where Zp„ and Z are given by 

Zpa = C^gpa = 2V(p^<^) and Z = Zp^g^" (60) 

(where () means symmetrization). 
Let us emphasize that the relation 

dpi^ = (61) 

is an identity which holds for all {gpv, (jpu, ^^). The vector field can be seen as 
generating an infinitesimal change of coordinates : contrarily to special relativity, one can 
build an infinite numbers of conserved vectors corresponding to all reparametrizations. 
These conserved vectors are useful to write integral conservation laws (known as "strong 
conservation laws"). When is a Killing vector field of the background, these strong 
conservation laws are called "Noether conservation laws", i.e. conservation laws related 
to a symmetry. 

We also stress that up to now, we have done no linearization; the case when is a 
"small perturbation" of Ai is considered in the next section. 
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5.2 Conserved vectors in a perturbed Friedmann-Lemaitre space- 
time 



We now apply this general formalism to perturbations in cosmology. The manifold is a 
perturbed Friedmann-Lemaitre universe and the background is taken to be a Friedmann- 
Lemaitre spacetime, as it is usual in the theory of cosmological perturbations. Let us note 
that other choices for the background can be made, for instance de Sitter spacetime (see 
e.g. |jl6|)- The line element of can then be written as 



ds'^ = g^^dx^dx" = {g^i, + h^u) dx^'dx'^ 



(62) 



where hnu is a "small perturbation'' 



We can now linearize equation (|56D . A first step yields 

1 



kIP 



5G>t + - {RW, 



R;sp]hp 



(63) 



The second step is to linearize Z^-. To do that we need to specify the arbitrary vector 
field As discussed before the relation 9^/'^ = is valid whatever field we use. However 
Killing vectors of the background (such that V(^^jy) = 0) are of special interest since they 
are related to symmetries and have a geometrical interpretation. Unfortunately the six 
Killing vectors of the Friedmann-Lemaitre background lay in spacelike hypersurfaces and 
cannot generate any quantity related to time translation, that is to energy. Now, we know 
that when two spacetimes are conformal, the Killing vectors of one of the spacetimes are at 
least conformal Killing vectors (that is such that ^{fj.^u) = jdiMu^ p^'^) of the other fl^. We 
also know that the Friedmann-Lemaitre spacetimes are conformal to Einstein, Minkowski 
and de Sitter spacetimes [0, each of which has a Killing vector field associated with 
time translation that induces a conformal Killing vector field on the Friedmann-Lemaitre 
spacetimes. In the following, we shall use one particular conformal Killing vector to 
generalize the symmetry under time translation. Let us first construct this vector field 
explicitely. 

The background metric being written as 



ds'^ = g^ydx^'dx" = a^{ri) —drf + •jijdx^dx^ with •jij = 5ij + K- 



^m.n.X X 



(64) 



we choose the conformal 



(65) 



where jij is the spatial metric in Weinberg's coordinates ||19 
Killing vector defined by 

e = ^o^ 

which is a Killing vector of Einstein static spacetimes, i.e. obtained by setting a = 1 in 
(p^). This vector is associated with translations in cosmic time, i.e. in the proper time 
of a comoving observer of the Friedmann-Lemaitre background. Other conformal Killing 
vectors could have been used (for instance the one associated with the conformity to de 
Sitter spacetime, see e.g. ^^), but (|65D is the useful one for the law we want to generalize. 

The vector having thus being chosen, we can now linearize Z^^. We see from equation 
(|60|) that Z^ vanishes for any Killing vector field. For any conformal Killing vector field 
equation (|59|) gives after linearization 0, 



SkZ^" = {hgpp - h^P) dpZ - ZV p {hgi^P - ^p) 
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with h = h^yQ^'^ and Z given by equation 

Let us now, as a third step, speciahze to scalar perturbations, that is to the sub-class 
of perturbations in synchronous gauges as defined in ([^). We now have all the elements to 
compute the conserved vector associated to the conformal Killing vector field (|65D. One 
can first easily show that 

and Z = Sn, (67) 



and a straightforward calculation then yields : 

4 



a 



n 



k6T^ + kQ^o + -^{D'hf - D^h) 



In synchronous gauge ( |68l) can be written in terms of E and h as 



7[«: {6Too + eoo) + Kh-nh'] 
K{6Tok + eok)+2n^''di(E 



-h- 



(68) 



(69) 



It is clear that the identity (0) with /'^ given by (p8| ) or (p9| ) is a consequence of the 
Einstein constraint equations (p!2|-p!3|), but a useful one as we shall see. 



5.3 



Encoding the matching conditions in a conserved vector 
field 



3), and recaUing that h = h~ 



2nr 



Comparing the expression of Too (p4 ) and tq^ ( |53D to 
on superhorizon scales, we have that, in that limit : 

- ^'l - 1 

"Too — OkT^k — 

Thus, for superhorizon scales, the matching condition (|3T| ) can be encoded in the compo 
nents of a conserved vector field related to a symmetry since 



(70) 



[t"00]: 







0. 



(71) 



This enables us to relate r^o^ in the long wavelength limit, to a vector field built with the 
conformal Killing vector field associated with time translation. In fact the construction 
of T^y uses in an obvious way the fact that a Friedmann-Lemaitre universe with fiat spatial 
sections is conformal to Minkowski spacetime (see Verraraghavan and Stebbins [|11[] section 
VI and Pen et al. ||] section IV-F). 

Too is often referred to as a "pseudo-energy" ||^. Our construction can justify such a 
denomination since it is related to a symmetry involving time translation. 
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Appendix 1 

We compare the terms c^A\l/ and (7i^ — K){1 + Scl)"^ in equation (|^) for modes of 
characteristic length L (i.e. such that A\I/ ~ in a universe with hyperbohc spatial 

sections {K = —1) and dominated by a mixture of radiation (subscript r) and dust 
(subscript m). 

Using the fact that P = Pj. and p = Pm + Pr, we can relate the sound speed = P'/ p' 
to the energy density of the radiation and the dust by 

4 Pr 

If we now introduce x = a/a^ and VLq = np^ /3Hq (where a subscript means that we 
evaluate the quantity today and X stands for either r or m), and use the Friedmann 
equation (P!7|), we have that 



clA^ « in' -K)il + 3cl)^ ^ « 3 f 1 + lx%\ (2+ 







(73) 



^ 4 ^5 ; V 1 - ^0 a;2 
All modes belonging to the categories (II) and (III) are such that 

L-'<R-hI. (74) 
Therefore, the inequality (|73D is satisfied for all x such that 



1 + ^ < 3 1 + -x-^ > ' ^ ' ° 



(75) 



For fi™ ~ 0.2, it is easy to see that this inequality is satisfied for all x in ]0, 1]. 
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